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We report results and analysis of time-resolved photoinduced reflectivity experiments on the 
cuprate superconductor YBa2Cu306.5. The sample, which has rc=45 K, was characterized by a 
high degree of purity and Ortho If ordering. The change in reflectivity Ai? was induced and probed 
using pulses of 100 femtosecond duration and photon energy 1.55 eV from a Ti:Sapphire laser. We 
provide a detailed picture of the decay rate 7 of A_R as a function of temperature T and pump 
intensity /. At low T, 7 decreases linearly with decreasing 7, extrapolating to nearly zero in the 
limit that / tends to zero. At higher temperature 7 has the same linear dependence, but with 
nonzero limit as / ^ 0. In the interpretation of these results we assume that A_R is proportional to 
the nonequilibrium quasiparticle density created by the laser. From an analysis of the 7 vs. / we 
estimate /3, the coefficient of proportionality relating the quasiparticle decay rate to the density. The 
intercept of 7 vs. / yields the thermal equilibrium quasiparticle decay rate. In a discussion section, 
we argue that the quasiparticles induced by the laser occupy primarily states near the antinodal 
regions of the Brillouin zone. We explain the divergence of the lifetime of these particles as T and 
I both tend to zero as a consequence of momentum and energy conservation in electron-electron 
scattering. Next, we discuss the significance of the measured value of /3, which is ~ 0.1 cm^s~^. We 
point out that the natural unit for /3 in a two-dimensional superconductor is h/m*, and define a 
dimensionless constant C such that /3 = Ch/m* . If the decay process is one in which quasiparticles 
return to the condensate with emission of a phonon, then C is a measure of the electron-phonon 
interaction. Alternatively, expressing the marginal Fermi liquid scattering in the normal state in 
terms of an effective /3 implies C = I/tt, which is in excellent agreement with the experimentally 
determined value in the superconducting state. 

PACS numbers: 74.25.Gz, 78.47.+P 



I. INTRODUCTION 

A special property of the cuprate superconductors is 
that the energy required for the creation of a quasipar- 
ticle depends on the direction of its momentum Jj. The 
creation energy is zero for momenta in the 'nodal' direc- 
tion, oriented at 45° relative to the Cu-0 bond. The most 
energetically expensive quasiparticles are the 'antinodal' 
ones, whose momenta are nearly parallel to the bond. 
Because they feel the pairing interaction most strongly, 
their properties may hold the key to high- Tc supercon- 
ductivity. Unfortunately, their tendency to form strong 
pairs makes them difficult to study. In thermal equilib- 
rium the population of quasiparticles is overwhelmingly 
dominated by the low energy nodal ones. As a result, 
transport measurements performed in equilibrium, such 
as microwave ^ and thermal Q conductivity, are insen- 
sitive to antinodal quasiparticles. 

A potentially powerful approach to studying the inter- 
actions of antinodal quasiparticles at low temperature is 
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to create a nonthermal population by external excitation. 
By probing the relaxation of the nonthermal population 
to the ground state, one may hope to learn about inter- 
actions between quasiparticles that would not normally 
be present at low temperature. Time-resolved optical 
techniques are ideally suited to creating a nonequilibrium 
density of quasiparticles and measuring their subsequent 
relaxation. These techniques are based on mode-locked 
lasers that produce pulses of 10-100 fs duration with a 
wide range of energy per pulse, center wavelength, and 
repetition rate. Such measurements are performed in a 
'pump-probe' mode in which a beam of pulses is split in 
two. The pump pulse creates the nonequilibrium state 
while the probe pulse senses the change in the optical re- 
sponse of the medium due to the nonequilibrium popula- 
tion. The time delay between the two pulses is controlled 
continuously and accurately by varying the optical path 
length difference between the two beams. Measuring the 
transmission or reflection of the probe as a function of 
time delay gives information about the return to equilib- 
rium after the pulsed excitation. 

There have been several time-resolved optical measure- 
ments performed on the YBa2Cu307_2: (YBCO) system 
of cuprate superconductors. The earliest work 0,|a,SI3 
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reported the change in the reflectivity (Ai?) of a 1.5 eV 
probe due to photoexcitation at the same energy. The 
measurements showed that in the normal state Ai? is 
small and decays very rapidly. Upon cooling below Tc, 
Ai? increases rapidly, and its decay rate decreases. These 
results suggested that carrier thermalization and/or re- 
combination proceed rapidly in the normal state, but are 
strongly impeded by the opening of the superconducting 
gap. 

Subsequent measurements provided a more detailed 
picture of the magnitude and decay rate of Ai? as a 
function of carrier concentration jsf . In underdoped sam- 
ples Ai? is readily detectable in the normal state. The 
decay rate slows with cooling, suggesting a correlation 
between the relaxation rate of the nonequilibrium state 
and the appearance of the pseudogap. It was reported 
that the relaxation time is insensitive to T in the su- 
perconducting state, with the exception of a peak near 
Tc- The results were interpreted in terms of an excited 
state in which quasiparticles and phonons rapidly reach 
quasiequilibrium. The decay rate was conjectured to re- 
flect the thermalization of nonequilibrium phonons. 

Recently, measurements of Ai? at 1.5 eV were reported 
on an untwinned single crystal of YBa2Cu306.5 Ortho II 
Ik]. (The des ignation Ortho II refers to the macroscopic 
ordering of the atomic layer containing the Cu-0 chains. 
In the Ortho II phase the chains alternate between fully 
occupied, and entirely unoccupied, by O atoms). In these 
measurements the range of pump intensity was extended 
to nearly two orders of magnitude smaller than used in 
earlier measurements. It was shown that while the T 
dependence of the decay rate is weak when the intensity, 
i, is large, it becomes extremely strong as I is reduced. 
In fact, the decay rate appeared to vanish as T and I 
both tend to zero. It was suggested that the strong I 
dependence of the decay rate is not consistent with a 
picture in which the quasiparticles and phonons reach 
quasiequilibrium and an alternate mechanism involving 
the pairwise scattering of quasiparticles was suggested. 

In this paper, we present further measurements and 
analysis of time-resolved photoinduced reflectivity in 
YBCO Ortho II. The experimental apparatus and sample 
characterization are described in Section II. In Section III 
we introduce the Rot hwarf- Taylor (RT) equations, which 
provide a phenomenological framework for interpreting 
nonequilibrium dynamics in superconductors. Section IV 
presents measurements of the decay of the photoinduced 
reflectivity as a function of time, temperature, and pump 
beam intensity. The essential experimental flnding, as in 
Ref. 9], is that the characteristic decay time diverges as 
both T and I approach zero. However, we present several 
results beyond those already reported. We present and 
analyze the time dependence of the transient reflectivity, 
showing that it is described well by the pairwise scatter- 
ing of quasiparticles. Second, we report improved mea- 
surements of the asymptotic value of the decay rate in 
the limit that / goes to zero. According to the RT equa- 
tions, the low- intensity limit is a direct measure of ^th^ 



the quasiparticle recombination rate in thermal equilib- 
rium. Through improvements in sensitivity we were able 
to measure ^th with more than one order of magnitude 
greater precision than previously. Finally, in Section V 
we present analysis and interpretation of the experimen- 
tal results. We argue that the quasparticles that give rise 
to Ai? are antinodal in character and interpret the decay 
rate as a measure of the strength of antinodal quasipar- 
ticle interactions. 



II. EXPERIMENTAL METHODS 

Pump and probe measurements were performed on a 
mechanically detwinned single crystal of YBa2Cu306.5 
with rc=45 K. The sample was grown in a BaZrOa cru- 
cible, which yields material with purity at the 0.99995 
level. The high purity allows allows the development of 
very long correlation lengths (jg = 148A, ^j, = 430A, 

= 58 A) of Ortho II order I1QI| , in which the the charge 
reservoir layer consists of alternating filled and empty 
copper oxygen chains. Together with YBa2Cu408, it is 
one of two underdoped cuprate superconductors in which 
doping does not introduce disorder. The relative lack 
of disorder is reflected in the low temperature transport 
properties, which indicate quasiparticle scattering times 
in excess of 30 ps [llj . 

In our experiments, the pump and probe beams were 
produced by a mode-locked Ti:Sapphire laser. The pulses 
have duration 100 fs, repetition rate 12 ns, and center 
wavelength 800 nm. Both pump and probe beams were 
focused onto the sample with a 20 cm focal length lens, 
yielding a spot size of 75 /xm diameter. The reflected 
probe beam was focused onto a Si photodiode detector. 

A double modulation scheme was used to optimize sen- 
sitivity to small changes in the power of the reflected 
probe beam. In order to minimize the noise from the Si 
photodiode detector, a photoelastic modulator was used 
to modulate the amplitude of the pump beam at 100 
KHz. In addition, a galvanometer-mounted mirror var- 
ied the path-length difference between pump and probe 
beams at a frequency of 40 Hz. This yields a rapid scan 
of the time delay between pump and probe, which helps 
to suppress noise due to 1/f fluctuations of the probe 
beam power. In order to demodulate and extract Ai? 
as a function of time delay, the output of the Si photo- 
diode was sent to a lock-in amplifler for phase-sensitive 
detection of the 100 KHz component. The output of 
the lock-in was then sent to a digital oscilloscope whose 
time base was triggered synchronously to the oscillating 
mirror. With the oscilloscope in averaging mode, a mini- 
mum detectable Ai?/i? of « 10~^ could be achieved after 
approximately 10 minutes of accumulation time. 

In this study, we focused on the decay rate of Ai? im- 
mediately following the pulsed excitation, or 7(0). For 
most of the measurements, we determined 7(0) by flt- 
ting Ai?(t) at small time delays by a decaying exponen- 
tial. However, special considerations arose in measure- 
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ments of 7(0) at very low pump power. As described 
in succeeding sections, 7(0) decreases as the pump in- 
tensity is lowered. At low temperatures, the decay rate 
decreases to the extent that the change in AR over the 
25 ps time delay range produced by the oscillating mirror 
is extremely small. To measure 7(0) in this regime, we 
switched to a detection scheme that probes directly the 
derivative of AR with respect to time delay. The signal 
from the first lock-in was sent to a second lock-in, rather 
than the digital oscilloscope, for phase-sensitive detec- 
tion at the frequency of the oscillating mirror. When the 
oscillating time delay is less than the decay time of AR, 
the output of the second lock-in is proportional to the 
derivative dAR/dt. In this mode of data acquisition we 
vary the pump-probe delay by a conventional system of 
a retroreflector mounted on a translation stage. 

When measured at low pump intensity, 7(0) is re- 
markably sensitive to T, decreasing about an order of 
magnitude between 15 and 10 K, for example. The ex- 
treme sensitivity to temperature suggests that laser heat- 
ing can introduce significant error in determining the 
T-dependence of 7(0). (By laser heating we refer to 
the steady-state increase in sample temperature due to 
the time-averaged laser power) . For the experiments de- 
scribed above, the pump intensity was lowered below the 
intensity of the probe. In this regime, the probe beam is 
responsible for heating the photoexcited region. At low 
temperature we can readily observe that 7(0) varies as 
the probe intensity, and consequently the sample tem- 
perature, is varied. The obvious remedy of reducing the 
probe power has the drawback that the AR signal soon 
disappears below the detector noise level. 

Empirically, we found that the optimal compromise be- 
tween signal-to-noise ratio and laser heating is obtained 
at a probe power of approximately 1.2 mW. At this power 
level 7(0) can be measured, yet the induced temperature 
change is significantly less than the bath temperature. 
Determining the dependence of 7(0) on T at this power 
level requires that we know the temperature of the pho- 
toexcited region of the sample. Unfortunately it is ex- 
tremely difficult to measure, with the required accuracy, 
the temperature of the sample directly under the focal 
spot. 

Because of this difficulty, we have performed a numer- 
ical simulation of the thermal diffusion equation to esti- 
mate the laser-induced temperature change. We input to 
this simulation a realistic model for the thermal proper- 
ties of the sample and its coupling to the bath. In the 
experiment, the sample is attached to a sapphire plate 
using thermal grease. After mounting, the thickness of 
the grease layer is determined using an optical micro- 
scope. The input parameters to the numerical simula- 
tion are the thicknesses of the sample and grease layer, 
as well as their temperature dependent thermal conduc- 
tivities 0, 01 . Figure 1 shows the temperature of the 
photoexcited region as a function of the bath tempera- 
ture, as determined by finite-element analysis of the dif- 
fusion equation. In the simulation, the laser power is 1.2 
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FIG. 1: Temperature of the photoexcited region of the sam- 
ple as a function of the bath temperature, as determined by 
numerical simulation. 



mW, focused to a spot of diameter 75 jjum. The thick- 
ness of the sample and thermal grease were 20 /zm and 5 
/im, respectively. In the subsequent analysis of the low- 
pump intensity data, we used the calculated temperature, 
rather than the temperature of the thermal bath. 



III. BACKGROUND: ROTH WARE- TAYLOR 
EQUATIONS AND THE 
PHONON-BOTTLENECK 

In this section we introduce the Rothwarf- Taylor 
equations, which provide a successful phenomenological 
framework for understanding nonequilibrium dynamics 
in superconductors Il5l |. The RT equations are a pair of 
rate equations that describe a system of superconduct- 
ing quasiparticles coupled to phonons. In the RT phe- 
nomenology the excited state is characterized by number 
densities rather than nonequilibrium energy distribution 
functions. This drastic simplification is justified in s-wave 
superconductors because quasiparticles rapidly thermal- 
ize to a narrow range of energy just above the gap. In 
d-wave superconductors the time evolution of the quasi- 
particle distribution function may be more complicated. 
In spite of this, we have found that the RT equations pro- 
vide an excellent description of the dynamics of YBCO 
Ortho II after pulsed photoexcitation. In the final section 
we comment on the underlying reasons for the applica- 
bility of the RT approach. 

In this work we write the RT equations in the following 
form. 



n = V -f 2Ar7pc - I3v? (1) 
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N - Iph + I3n^l2 - jpcN - (iV - Neq)iesc (2) 

where n, and A^, are the number densities of gap en- 
ergy quasiparticles and phonons, respectively. The right 
hand side of Eq. 1 expresses the difference between the 
rates of quasiparticle creation and annihilation. Iqp is the 
external generation rate and 2N"fpc is the rate of pair 
creation via annihilation of gap energy phonons. The 
quasiparticle annihilation rate /3n^ varies quadratically 
with density because recombination is a two-particle scat- 
tering event. In Eq. 2 the time rate of change of the 
phonon density is given by the external phonon creation 
rate and the same recombination and pair creation terms 
(with opposite signs) as in Eq. 1. The parameter 7esc, 
which appears in the last term in Eq. 2, is the rate 
at which a gap energy phonon decays into below-gap en- 
ergy phonons that cannot regenerate a quasiparticle pair. 
Typically, j^sc ^ 7pc, in which case a gap energy phonon 
is far more likely to regenerate a quasiparticle pair than 
to decay into the bath. Although the escape process is 
relatively slow, it is essential for the ultimate return of 
the coupled system to equilibrium. 

In the present paper, we are concerned with the evolu- 
tion of the system following pulsed excitation, in which 
case there are two characteristic time regimes. For de- 
lay times t such that t ^ 7pc^, the quasiparticles and 
phonons have not yet reached quasiequilibrium. In this 
regime the quasiparticle population decays at the density 
dependent rate /3n. For t ^ the exchange of energy 
between the phonons and quasiparticles has brought the 
two populations into quasiequilibrium. In this regime N 
and n both decay at the much slower rate 7esc- The de- 
cay of the quasiparticle population has become slow and 
independent of density, despite the fact that the rate of 
quasiparticle scattering may be rapid and density depen- 
dent. The limit on the decay of the nonequilibrium quasi- 
particle density imposed by the quasiequilibrium with the 
phonons is termed the phonon bottleneck. 



IV. EXPERIMENTAL RESULTS 

A. Intensity dependence at low temperature 

The decay rate of the photoinduced reflectivity in 
YBCO Ortho II depends strongly on both the temper- 
ature, T, and the pump intensity, /. Fig. 1 illustrates 
the decrease of the decay rate with decreasing / at a 
fixed T=9 K. Each curve, measured using a different / 
in the range from 2 to 25 mW, has been normalized to 
the same value at t = to illustrate the variation in de- 
cay rate. Although the decay rate continues to decrease 
as / is lowered below 2 mW, these curves are not shown 
because their noise fluctuations obscure the higher inten- 
sity curves. We discuss the decay rate at very low pump 
intensity in more detail in Section IV. C. 
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FIG. 2: AR/R at T=9 K as a function of time following- 
pulse excitation, for several pump intensities. All curves were 
normalized to the same value at delay time zero to illustrate 
the variation in initial decay rate. The most rapid decay is 
seen at the highest pump intensity, corresponding to 25 mW 
of laser power focused to spot of diameter 75/^m. The decay 
becomes systematically slower as the intensity is reduced by 
factors 0.69, 0.44, 0.24, 0.15, 0.10, and 0.06. 



The variation in decay rate with excitation density ob- 
served in YBCO Ortho II is highly unusual. Analogous 
experiments on s-wave superconductors have consistently 
found that the decay rate of the nonequilibrium state is 
independent of the quasiparticle density [l5j. The lack 
of density dependence is understood as a manifestation 
of the phonon bottleneck effect. As discussed in the pre- 
vious section, the bottleneck sets in when t reaches ap- 
proximately . If there is no significant recombination 
on this time scale, that is if /3n ^ 7pc, then only density- 
independent bottleneck dynamics can be seen. 

Within this picture, the absence of a bottleneck in 
YBCO Ortho II indicates that the hierarchy of pair cre- 
ation and recombination rates is reversed as compared 
with conventional superconductors. (In Section V.D we 
discuss why the hierarchy of rates may be reversed in 
cuprate superconductors). If /3n ^ jpc, then in the time 
window t < jp^ quasiparticle recombination proceeds 
in the absence of the reverse process of pair creation. 
Measuring the decay of the quasiparticle density in this 
time window provides direct information about the rate 
of quasiparticle-quasiparticle scattering. 

To test whether our measurements of YBCO Ortho 
II are in a pre-bottleneck regime, we have analyzed the 
time-dependence of the decay of AR. Fig. 3 is a double 
log plot of AR as a function of time delay, for several 
pump intensities. At the highest excitation density AR 
decays as a power law for t > 3 ps. With decreasing 
excitation density the onset of power law decay shifts to 
longer times. At the lowest intensity the crossover time 
is greater than the maximum time delay in this set of 
measurements, which is ^25 ps. 
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Eq. 5 predicts that the excited population approaches 
1/Pthit l/(3nph{0), regardless of the initial quasipar- 
ticle density. However, as is clear from Fig. 3, the mea- 
sured decay curves do not obey this prediction. At long 
times the decay is closer to t~^-^ than t~^, and curves 
for different intensities do not merge to a single curve at 
long times. 

The discrepancies described above can be traced to a 
complication omitted from the preceding analysis. Be- 
cause the intensity of the pump beam decreases expo- 
nentially with increasing depth below the sample surface, 
z, the local quasiparticle density, nph{z,t), is spatially 
nonuniform. The initial local density np/i(z,0), equals 
f^p/i(0, 0)e~"^, where a is the absorption coefficient at 
the pump wavelength. Assuming negligible diffusion of 
excitations in the z direction, the local excitation density 
decays as, 



nph{z,t) = 



nph{z,0) 



[1 + I3nph{z,0)t] 



(6) 



FIG. 3: AR/R vs. time delay at T=9K plotted on a double 
logarithmic scale, for pump laser powers (in mW) 0.1, 0.2, 
0.4, 0.6, 1.0, 1.5, 2.4, 3.8, 6.0. 



We compare these curves with the predictions of the 
RT equations in the time regime before the onset of the 
phonon bottleneck. This decoupled regime corresponds 
to the 7pc limit. In this limit the form of the RT 
equations is clarified if the total quasiparticle population 
is written as Uph + nth, where Uph and nth are the pho- 
toinduced and thermal equilibrium quasiparticle density, 
respectively. Similarly, we define Nph and Nth as the 
photoinduced and thermal equilibrium phonon popula- 
tions. Detailed balance relates the thermal densities of 
phonons and quasiparticles such that, 



The probe beam, whose change in reflectivity measures 
the excited quasiparticle density, also decays exponen- 
tially in the sample. Therefore the measured reflectivity 
change is related to an exponentially weighted average of 
the nonequilibrium quasiparticle density p^ . 



M{t)=a-^ dze-"'nph{z,0)/[l + l3nph{z,Qi)t]- (7) 



Assuming that the change in reflectivity is proportional 
to N finally yields. 



Ai?(t) 



2Ai?(0) ln{l+-iQt) 



lot 



-[1 



lot 



(8) 



'^Nthlpc = Pnl^. 



(3) 



Substituting the above relation into Eq. 1, and setting 
Ipc — 0, yields a rate equation for Uph that is independent 
of the phonon population, 



"^ph — Iqp P^ph 



'2(3nthnph- 



(4) 



Eq. 4 is equivalent to the rate equation applicable to 
nonequilibrium electrons and holes in semiconductors, 
and has been studied extensively in connection with pho- 
toconductivity 13}- If nth vanishes as T approaches zero, 
then the rate equation approaches fiph = Iph — l^n^h ■ 
low T the quasiparticles obey simple second order, or bi- 
molecular, reaction kinetics. Integration yields the decay 
after pulsed excitation, 



riphit) 



nph{0) 



[1 + I3nph{0)t] 



(5) 



where 70 = PnphiO, 0). 

The solid lines in Fig. 3 indicate the best fit of Eq. 8 
to the experimental data. The fit is obtained by varying 
Ai?(0) and 70 for each curve. There is excellent agree- 
ment between the data and the prediction of the bimolec- 
ular rate equation when the exponential variation of the 
pump and probe intensities is taken into account. The 
agreement suggests that for t < 25 ps, the phonon bot- 
tleneck is not yet established and the decay rate of AR is 
a direct measure of the quasiparticle-quasiparticle scat- 
tering rate. 

From the fitting procedure we determine Ai?(0) and 
70 for each curve. In Fig. 4 we plot 70 as a function 
of AR{0)/R. It is evident that the characteristic decay 
rate increases linearly with the magnitude of the initial 
reflectivity change. The linear relationship indicates that 
the entire family of decay curves can be described by a 
single bimolecular coefficient (3. 
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FIG. 4: The initial decay rate at the sample surface, 70, as 
a function of the initial reflectivity change, AR{0)/R. The 
dashed line emphasizes the linear dependence of initial decay 
rate on initial density. 



B. Temperature dependence 



FIG. 5: AR/R induced by 2.5 mW of laser power, as a func- 
tion of time following pulse excitation, for several tempera- 
tures in the range from 5-70 K. All curves were normalized to 
the same value at delay time zero to illustrate the variation 
in initial decay rate. The decay is most rapid at the highest 
temperature, and becomes systematically slower as T is re- 
duced from 70 K to 50 K, 44 K, 28 K, 22 K, 17 K, 12 K, and 
5 K. 



In the previous section we showed that the decay of 
AR at low temperature can be described by the bimolec- 
ular rate equation with a single value of the scattering 
coefficient (3. In this section we analyze the effect of rais- 
ing the temperature on the rate of decay of AR. Fig. 
•5 shows a scit of normalized decay curves of AR, mea- 
sured in the temperature range from 5-70 K, induced by 
a pump power of 2.5 mW. The set of curves illustrates 
clearly that the decay rate at fixed / increases rapidly 
with increasing T. 

To analyze the T dependence, we again consider the 
RT (equation in the decoupled regime (Eq. 4) . At nonzero 
T we must include the last term on the right-hand side, 
which was neglected previously. Physically, this term 
describes the rate of scattering of a photoinjected quasi- 
particle by a thermal equilibrium one. 

We have found that the most direct way to compare 
the data with Eq. 4 is to focus on the initial decay rate, 
7(0), of the experimentally measured transients. Here 
7(0) is defined as the limit of the instantaneous decay 
rate ■y{t) = —hph/nph as t approaches zero. According 
to Eq. 4, 7(0) — /3[np/i(0) + 2nth\- Assuming that the 
initial reflectivity change is proportional to np/i(0), a plot 
of 7(0) vs. Ai?(0) should yield a straight line with slope 
proportional to /3 and intercept 2(5nth- 

Fig. 6 shows 7(0) plotted as a function of Ai?(0) / R for 
several temperatures in the superconducting state. The 
data at the lowest temperature are essentially equivalent 
to the results presented earlier. In Fig. 4 the decay rate 
was determined by fitting the entire time dependence, 
while in Fig. 6 the decay rate is determined from the 
initial slope. The reason the decay rates in the two figures 
differ by a factor of two can be understood from Eq. 8. 



Taking the derivative of this formula with respect to time 
shows that 7(0) = 7o/2. 

With increasing T the 7(0) vs. Ai?(0) plots shift verti- 
cally, with little change in slope. The fact that the slope 
is nearly constant implies that (3 depends weakly, if at 
all, on the temperature. The increase of the intercept 
with T implies a rapidly increasing density of thermal 
equilibrium excitations. In the next section we examine 
the temperature dependence of the intercepts, as deter- 
mined from measurements performed at very low pump 
and probe intensities. 

Fig. 7 presents a different perspective of the decay 
rate as a function of temperature and intensity: 7(0) 
as a function of T for three values of /. In the nor- 
mal state and just below Tc the different data sets lie 
on the same curve, indicating that 7(0) is independent 
of /. As the temperature is lowered below Tc the decay 
rates measured with different pump intensities begin to 
diverge. As the temperature tends to zero, the decay 
rate crosses over to a temperature independent regime, 
with the crossover temperature higher when the rates are 
measured at greater pump intensity. A consistent expla- 
nation for the T and / dependence of 7(0) can be found in 
the picture of pairwise scattering involving both thermal 
equilibrium and photoinduced quasiparticles. The lack 
of intensity dependence above Tc suggests that in this 
regime the thermal density of quasiparticles is far larger 
than the photoinduced density. The strong temperature 
dependence above Tc suggests that nth decreases with 
decreasing temperature in the normal state, possibly due 
to the opening of the pseudogap. The onset of inten- 
sity dependence at Tc may indicate that nth has become 
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FIG. 6: Initial decay rate, 7(0), as a function of initial reflec- 
tivity change, AJ?(0)/iZ, for T=12 K, 17 K, 22 K, and 30 K. 
The thermal equilibrium decay rate, obtained from extrap- 
olation of the data to zero A_R(0)/-R increases rapidly with 
temperature. 
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C. Low intensity regime 

In this section we focus on measurements of 7(0) per- 
formed at the lowest laser intensities that are accessible 
experimentally. The motivation for studying the low- 
intensity regime is to measure, via a nonequilibrium ex- 
periment, the recombination lifetime of quasiparticles in 
equilibrium. In order to probe equilibrium properties, 
the density perturbation introduced by the laser must be 
small, such that Uph ^ nth- According to Eq. 4, 7(0) 
approaches 2j3nth in this limit, which is exactly twice the 
thermal equilibrium recombination rate, ^th- 

Measuring ^th accurately becomes increasingly diffi- 
cult at low temperature because of the requirement that 
riph ^ nth- As nth decreases rapidly with decreasing T, 
the pump intensity must be lowered accordingly. Eventu- 
ally Ai? reaches the noise floor of the experiment, which 
sets a limit on the smallest ^th that can be measured. 
In principle, the signal size could be increased by rais- 
ing the intensity of the probe beam. However, this in- 
evitably leads to an increase in the average temperature 
of the photoexcited region compared with the bath tem- 
perature. This problem becomes increasingly severe at 
low temperature, where the thermal conductivity of the 
sample is small. For the data to be presented, the probe 
intensity was maintained at 1.2 mW, a value which pro- 
vides the best compromise between sensitivity and heat- 
ing. 

At each T, we obtain 274/1 from the extrapolation to 
zero of a linear fit to 7(0) vs. Ai?. The temperature 
of the photoexcited region is determined through the nu- 
merical analysis described in Section II. Fig. 8 is a double 
logarithmic plot of 274/1 as a function of T, showing re- 
sults for both YBCO Ortho II and a thin film sample 
of BiaSrsCaCuzOg (BSCCO) [H). At the upper limit of 
the T range, ^th is nearly identical in the two materials. 
However, as mentioned previously, the dependence of 'yth 
on T becomes very strong in YBCO Ortho II at low tem- 
perature. If the functional form was assumed to be a 
power law, then ^th would approach a dependence. 

We believe it to be more reasonable to assume that 
^th in YBCO Ortho II decreases exponentially at low 
temperature. To analyze this exponential dependence, 
we compare the data with the formula. 



FIG. 7: Initial decay rate, 7(0), as a function of T for three 
values of pump power. 



comparable to nph- Alternatively, the sudden change in 
decay kinetics may be related to the onset of the coher- 
ence in the antinodal quasiparticle self-energy. Finally, 
the decay rate measured at fixed pump intensity crosses 
over to T-independence when nth becomes much smaller 
than nph{T). The overall behavior of 7(0) suggests that 
the decay rate vanishes in the limit that both / and T go 
to zero. 



2lth{T) = -^—exp[~Ath{T)/kBT] . (9) 
n Jo 

We introduce a prefactor proportional to for two rea- 
sons, one theoretical and the other experimental. The 
theoretical reason is that nth{T) is predicted to vary as 
for a d-wave superconductor. The experimental rea- 
son is that in some samples of BSCCO 74/1 is propor- 
tional to over the temperature range from 10 K to 
40 K [3. (The dashed hne in Fig. 8 indicates a 
dependence). Fitting Eq. 9 to the YBCO Ortho II data 
yields To = 1100 K and a temperature dependent ac- 
tivation energy A4;i that is plotted in the inset. The 
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FIG. 8: Twice the thermal equihbrium recombination rate, as 
determined by the extrapolation of 7(0) to zero laser power, 
as a function of temperature. Solid squares are results for the 
YBCO Ortho II sample, open circles show results obtained 
on a thin film sample of BSCCO with T^^ 85 K. The dashed 
line indicates a temperature dependence proportional to . 
The inset shows the activation energy of the recombination 
rate in YBCO Ortho II as determined from a fit to Eq. 9. 



origin of the exponential cutoff in ^th that appears in 
YBCO Ortho II and not in BSCCO is not understood 
at present. However, we note that analogous behavior is 
observed in the scattering rates observed by ac conduc- 
tivity experiments. In the YBCO system, the mean free 
path increases exponentially with temperature towards 
a limit of several microns set by elastic processes [Tl| . 
In BSCCO the T dependence is much weaker and the 
limiting value is only a few hundred A |l9j| . It is possi- 
ble that momentum conservation imposes constraints on 
the rate of quasiparticle-quasiparticle scattering in the 
cleaner YBCO system. These constraints may dictate 
that only Umklapp processes are allowed at low T, lead- 
ing to an exponential T dependence [23, 0| ■ 



V. DISCUSSION 

In analyzing the experimental data presented above, 
we begin (in Section V.A) with the most basic question: 
what type of excitation is probed by transient reflectiv- 
ity experiments performed on the cuprate superconduc- 
tors? We argue that these excitations are quasiparticles 
that occupy states close to the antinodal regions of the 
Brillouin zone. In Sections V.B and V.C we turn to a 
quantitative analysis of the decay rate of the photoex- 
cited state. Section V.B focuses on our observation that 
the lifetime of a photoinjected quasiparticle diverges in 
the limit that T and / both tend to zero. We show the- 
oretically that the constraints of momentum and energy 
conservation prevent thermalization of antinodal quasi- 
particles toward the nodes. In Section V.C we estimate 



the magnitude of the recombination coefficient [3 from the 
experimental data, and discuss its theoretical interpreta- 
tion. Finally, in Section V.D, we comment on possible 
reasons for the absence of a phonon bottleneck in the 
decay of the photoexcited state. 



A. What types of excitations are probed? 

As a first guess, one might imagine that the pho- 
toinduced quasiparticles occupy states near the nodes, 
as these are the lowest energy excitations. This intu- 
ition is consistent with the generally accepted picture of 
the quasiparticle distribution function in the nonequilib- 
rium state. In this picture the nonequilibrium quasipar- 
ticles rapidly adopt a Fermi-Dirac distribution, but with 
a chemical potential, and temperature, fi* , and T*, re- 
spectively, that are larger than their values in thermal 
equilibrium [2^ |2^ . In a d-wave superconductor this 
distribution describes a population of quasiparticles that 
is dominated by states near the nodes ('nodal' quasiparti- 
cles). The occupation of states with energy greater than 
max(/i*,T*) is exponentially small. 

On the other hand, the experimental evidence suggests 
that the photoexcited state is not a degenerate gas of 
nodal quasiparticles. Three observations, in particular, 
lead us to this conclusion: 

(1) The diffusion coefficient of photoinjected quasipar- 
ticles is at least two orders of magnitude smaller than 
the diffusion coefficient of nodal quasiparticles p3 | . The 
diffusive propagation of nonequilibrium quasiparticles in 
YBCO Ortho II was recently measured by the transient 
grating technique. Values for the diffusion coefficient of 
20 cm^/s and 24 cm^/s were determined for motion along 
the a and b crystalline axes, respectively. This may be 
compared with a nodal quasiparticle diffusion coefficient 
of « 6000 cm^/s as determined from microwave spec- 
troscopy on the same sample 'll] . 

(2) The reduction in the condensate spectral weight de- 
pends linearly on the pump intensity. Information about 
the photoinduced reduction of condensate spectral weight 
comes from visible pump-terahertz probe experiments 
[2^12^. At frequencies ~1 THz, the optical conductiv- 
ity 17(0;) is dominated by its imaginary part, a2, which is 
proportional to the condensate density. Following pho- 
toexcitation at 1.5 eV, 172 drops rapidly as the condensate 
density is diminished. At intensities below the saturation 
level, the loss of condensate density is a linear function 
of the energy deposited in the photoexcited volume. 

The linear reduction of condensate density with pump 
energy is inconsistent with the physics of nodal quasi- 
particles. This point can be illustrated by considering 
the temperature dependence of the superfluid density, 
Ps- Nodal quasiparticles are responsible for the linear 
in T decrease of ps observed in clean cuprate supercon- 
ductors [23. However, the energy stored in the nodal 
quasiparticle gas increases not as T, but as T^. There- 
fore the reduction in pa due to nodal quasiparticles varies 
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as the one-third power of the total energy in the quasi- 
particle gas. We would expect such a strongly sublinear 
relationship for a degenerate gas of nonequilibrium nodal 
quasiparticles as well. 

(3) Far-IR measurements of the photoexcited state 
show that the optical spectral weight that is removed from 
the condensate by photoexcitation shifts to energies on 
the order of 50-100 meV, or several times the value of 
the maximum gap p^ . As discussed above, thermal ex- 
citation of nodal quasiparticles does indeed remove spec- 
tral weight from the condensate. However, the spectral 
weight shifts to a very narrow (<30/j,eV) Drude peak cen- 
tered on w = 0. In contrast, the spectral weight removed 
from the condensate by photoexcitation shifts to frequen- 
cies that are wlOOO times higher than the width of the 
nodal quasiparticle Drude peak. This is further evidence 
against a degenerate gas of nodal quasiparticles. 

On the basis of the preceding arguments, we conclude 
that nodal quasiparticles do not dominate the population 
of quasiparticles created by photoexcitation. In other 
words, the quasiparticle distribution in the photoexcited 
state cannot be described as a Fermi-Dirac function with 
effective parameters ^* and T* . The experiments are 
more consistent with a distibution function that is peaked 
at an energy above the chemical potential. Such a dis- 
tribution can arise if it is not possible to scatter into the 
the nodal states during the lifetime of the photoexcited 
state. In Section 5.2 we show that the constraints of 
energy and momentum conservation severely restrict the 
rate at which hot quasiparticles scatter into states near 
the nodes. 



B. Stability of an isolated photoexcited particle 

In section V.A, we argued that the quasiparticle dis- 
tribution function in the photoexcited state peaks at an 
energy above the chemical potential and therefore is not 
of the Fermi-Dirac form. Such a nonequilibrium distribu- 
tion can arise if quasiparticles cannot scatter into nodal 
states during the lifetime of the photoexcited state. How- 
ever, we have also shown that this lifetime can be very 
long, reaching 0.6 ns at the lowest temperature and 
intensity accessible experimentally. By comparison a 30 
meV excitation in a Fermi liquid (with Fermi energy ^ 
1 eV) would decay into electron-hole pairs in ^3 ps. In 
the following we describe how the relative stability of a 
quasiparticle in a d-wave superconductor can arise from 
the severely restricted phase space for decay. 

In a Fermi liquid, quasiparticles readily decay into 
particle-hole pairs because energy and momentum can 
be conserved in the process. However, the kinematic 
contraints are much more severe in the case of a d- 
wave superconductor. Consider the momentum-resolved 
particle-hole excitation spectrum in the two cases. In the 
Fermi liquid the particle-hole spectrum forms a broad 
continuum that extends from zero wavevector to twice 
kp, even as the pair energy tends to zero. 



The d-wave superconductor differs from the Fermi liq- 
uid in that the Fermi contour shrinks to four nodal points. 
As a consequence, the pair excitation spectrum is much 
more localized in momentum space. Fig. 9 illustrates 
the momentum-resolved pair excitation spectrum of a d- 
wave superconductor. Each plot is a color-scale depiction 
of the density of two-particle excitations of a given energy 
in the first Brillouin zone. The spectra were calculated 
using a parameterization of the quasiparticle dispersion 
obtained by Norman 29] from a fit to ARPES data. At 
very low energies, the pair spectrum is localized near nine 
discrete wavevectors that represent the creation of par- 
ticles at the nodes. As the energy of the pair increases, 
the locus of possible wavevectors spreads out. However, 
for energies that are less than Aq, the allowed momenta 
remain clustered near the nine zero-energy wavevectors. 

In order for a single quasiparticle to scatter and emit 
a pair, its change in energy and momentum, Ae and Afc, 
must match that of the pair, epair and kpair- The con- 
straints imposed by this condition are illustrated in the 
diagrams shown in Fig. 10, which describe the possible 
decay events of a hot quasiparticle. For this example 
we have chosen a quasiparticle with energy 40 meV and 
momentum on the Fermi contour of the normal state. 
Each diagram corresponds to a fixed energy transfer in 
the scattering event. The color scale indicates the lo- 
cation of allowed values of Afc and kpair at the transfer 
energy. A scattering event is kinematically allowed only 
at momentum transfers where the two sets of wavevectors 
overlap. 

Fig. 10 illustrates how kinematics prevent the inelas- 
tic scattering of the hot quasiparticle. For example, con- 
sider the scattering events with the largest possible en- 
ergy change, 40 meV. To give up all its energy the hot 
quasiparticle must scatter to one of the nodes, thus there 
are only four allowed values of Afc. (The four points are 
not visible in Fig. 10(a)). One of these four wavevectors 
just touches the perimeter of the elliptical region of pair 
wavevectors, which indicates an allowed scattering event. 
In this scattering event the hot quasiparticle scatters to 
the node, with the simultaneous creation of a quasihole 
at the node and a quasiparticle at the energy and mo- 
mentum of the hot quasiparticle. This exchange process, 
although kinematically allowed, does not lead to a change 
in the particle distribution function. 

If the hot quasiparticle does not lose all of its energy, 
there are more choices for the final momentum. The locus 
of allowed Afc expands, becoming visible as arcs in panels 
(b) through (e) of Fig. 10. However, as the arcs expand 
the locus of pair creation wavevectors shrinks. Fig. 10 
shows that the region of allowed kpair shrinks faster than 
the expansion of allowed Afc, so that the two regions fail 
to overlap at any energy transfer less than the energy 
of the hot quasiparticle. The absence of overlap demon- 
strates that the thermalization of the hot quasiparticle 
is forbidden because of kinematic constraints. Recently, 
Howell et al. [IJ reached a similar conclusion based on 
an analytical, rather than numerical, approach. 



10 





(c) 






■ 


\ 


■ 




m 


\ 


m 


/ 


(f) 



FIG. 9: Color scale depiction of quasiparticle pair momenta in the first Brillouin zone for different values of the pair energy; 
(a) 10 meV, (b) 20 meV, (c) 30 meV, (d) 40 meV, (e) 50 meV, and (f) 60 meV. 



The kinematic constraints in this example are specific 
to quasiparticles near the normal state Fermi contour. 
They are a consequence of the rapid decrease of quasi- 
particle velocity with increasing distance from the node. 
For this reason, we anticipate that not all quasiparticles 
of a given energy will be stable. The contours of con- 
stant quasiparticle energy are distorted ellipses centered 
on each of the nodes. Quasiparticles near the major axis 
vertex of the ellipse will be stable, as we have seen. How- 
ever, the velocity for particles near the minor axis vertex 
is comparable the nodal quasiparticle velocity. For these 
quasiparticles there are no kinematic constraints to pre- 
vent rapid thermalization. Thus, we anticipate that after 
a short time following pulsed injection only nonequilib- 
rium particles near the major axis vertex will survive. 



C. Bimolecular recombination 

1. Estimating /3 

In the previous section we have seen that kinematic 
constraints can stabilize an isolated antinodal quasipar- 
ticle. However, our experiments show that a pair of 
quasiparticles can convert to a state which no longer con- 
tributes to AR. We can be confident that this conversion 
requires two particles, and not more, because the decay 
rate increases linearly with density. The coefficient that 
relates the decay rate to the density, /3, is directly related 
to the cross-section for inelastic scattering of two quasi- 
particles. Thus the magnitude of f3 is of fundamental 
importance as a measure of the coupling of quasiparti- 
cles to some other excitation of the interacting electronic 
system. 

We turn next to estimating /3 from the measured de- 
pendence of the decay rate on density. Experimentally, 
we measure both the nonequilibrium decay rate (3nph 
and the thermal equilibrium rate (3nth- To extract /3 
from these measurements, we need to know either Uph 
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FIG. 10: Depiction of the lack of overlap between the scattering momentum Afc and the pair momentum kpair- For this 
illustration the hot quasiparticle starts with energy 40 meV and momentum on the Fermi contour of the normal state. The 
arc segments in each square show the possible values of Ak for several values of the scattering energy Ae: (a) 40 meV; (b) 30 
meV; (c) 20 meV; (d) 10 meV; (e) 5 meV. The color scale plot in each of the squares shows the allowed values of kpair at the 
same energy. As the energy transfer decreases and the arc segments expand, the regions of kpair shrink faster. The lack of 
overlap is a graphical demonstration of how momentum and energy conservation prevent thermalization of a quasiparticle near 
the antinodal region of the Brillouin zone. 



or nth- Unfortunately, neither quantity is directly deter- 
mined from the data, and we need additional assumptions 
to estimate these quantities. To estimate Uph we use en- 
ergy conservation to convert from laser energy density to 
quasiparticle density. In the linear regime of 7 vs. / a 
recombination rate of 10^^ results from a laser energy 
density of 0.1 J-cm~'^. If photon energy is converted en- 
tirely to 40 meV quasiparticles, then their initial density 
is « 10""^^ cm~^, and therefore /3 ~ 0.1 cm^/s. 

We can also determine p from the thermal equilibrium 
recombination rate as well, if we know nth{T). For a 
d-wave superconductor, 

^ S (If )^ <-) 

where vp and V2 are the nodal quasiparticle velocities 



perpendicular and parallel to the Fermi contour, respec- 
tively 30]. Based on the RT equations, we would expect 
the thermal equilibrium decay rate to be proportional to 
TLth and thus vary with temperature as T^. Instead, we 
find experimentally that in YBCO Ortho II the limit- 
ing behavior of ^th is exponential rather than power law. 
The unexpected exponential cutoff of the recombination 
rate may be related to a similar exponential cutoff seen 
in the momentum scattering rate as measured through 
microwave spectroscopy 0. 

Despite the exponential cut-off of the recombination 
rate at low T, we can attempt to estimate (3 from the 
thermal recombination rate at higher T. Above k, 20 
K, jth{T) is not inconsistent with a quadratic depen- 
dence on temperature. In that range we found that 
'2lth{T) « {kBT/h)T/To, with To =1100 K. Assuming 
that 7t?i = (3nth, and using Eq. 10 to estimate nth(T) 
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(with vfIv-2, = 10) we find /3 = 0.26 cm^/s. Considering 
the approximations and assumptions involved, this value 
is quite close to the one determined from the nonequilib- 
rium recombination rate. 



2. Theoretical approaches for /3 

Within the context of BCS theory, the process of re- 
combination is the annihiliation of a quasiparticle pair 
with the simultaneous emission of a gap energy phonon. 
The theory of the recombination rate of BCS quasiparti- 
cles was developed by Kaplan et al. [sj, whose calcula- 
tions assumed an isotropic s-wave gap and the breaking 
of momentum conservation through strong elastic scat- 
tering ('dirty limit'). Their result for the recombination 
rate, expressed through our definition of recombination 
coefficient, is (3 = 4:7ra^{2A)F{2A)/hN{0), where a, F, 
and N, are the electron-phonon coupling strength, den- 
sity of phonon states, and density of electron states, re- 
spectively. While it is very useful to have this explicit 
result to compare with, we must bear in mind that the 
cuprate superconductors do not satisfy the assumptions 
upon which this calculation is based. The cuprates are 
2D rather than 3D, d-wave rather than s-wave, and are in 
the clean, rather than the dirty, limit of elastic scattering. 

We are not aware of a calculation of the phonon- 
mediated recombination rate applicable to the cuprate 
superconductors. However, in anticipation of such a 
calculation, we define a dimensionless constant C = 
hPN{0), that parametrizes the electron-phonon coupling 
strength in the clean limit. If in addition we express 
the 2D electronic density of states through an effective 
mass m*, then we obtain the simple formula (3 = CH/m*. 
Assuming an effective mass that is three times the free 
electron mass, h/m* = 0.3 cm^/s, so that (3 = 0.1 cm^/s 
corresponds to C w 1/3. It is striking that C is of or- 
der unity. This suggests that, while we do not know 
if phonon emission is indeed the mechanism for quasi- 
particle recombination, we have identified h/m* as the 
'natural unit' for (3. 

The body of experimental results on the cuprate su- 
perconductors have led many researchers to conclude 
that quasiparticle scattering is determined by electron- 
electron, rather than electron-phonon interactions. The 
BCS theory of recombination was adapted by Quinlan 
et al. |33| to the case where quasiparticle-pair energy is 
converted to antiferromagnetic spin fluctuations, rather 
than phonons. As we explain below, the results of this 
calculation cannot be compared directly with /?. The re- 
combination rate at low T is predicted to vary as T^. 
Therefore 74/1 is not proportional to nth, which, as we 
have seen, varies as T^. Expressed in the language of 
this paper, /3 vanishes in the limit that T goes to zero, 
in contrast to the T-independent value observed experi- 
mentally. 

The prediction that 7 cx T'^ while nth oc is a re- 
flection of the underlying Fermi liquid theory (FLT). In 



FLT the normal state scattering rate varies as T^, while 
the thermal quasiparticle density varies as T. The cor- 
responding quantities vary as one higher power of T for 
the d-wave superconductor because its density of states is 
linear in energy, rather than contant. Thus it is a general 
consequence of the phase space of Fermi liquids that the 
quasiparticle scattering rate is not simply proportional 
the quasiparticle density. 

It is interesting, and perhaps relevant, to recall here 
the celebrated fact that the high-T^ cuprates violate this 
central prediction of FLT. That is to say, the normal 
state scattering rate varies linearly, rather than quadrati- 
cally, with the temperature. In this marginal Fermi liquid 
(MFL) 33], the scattering rate ~ 2kBT/h is proportional 
to the quasiparticle density N{0)kBT. The coefficient of 
proportionality, or Pmfl, is simply h/nm*. Evaluating 
this expression with m* — SniQ yields cm^/s, 
which is remarkably close to the experimental value. This 
agreement suggests that there may be a close connection 
between the quasiparticle scattering in the normal state 
and the antinodal quasiparticle recombination in the su- 
perconducting state. 

D. Absence of phonon bottleneck 

In the previous section we described two mechanisms 
for quasiparticle pair recombination. In the more famil- 
iar one, the pair is destroyed with simultaneous creation 
of a 2A phonon. In the second mechanism the pair is 
converted to an electronic excitation, such as an antifer- 
romagnetic fluctuation. Both mechanisms can provide 
an explanation for the absence of a bottleneck. In the 
second mechanism, the quasiparticle pair decays to other 
particle-hole excitations. If the number of particles gen- 
erated in the pair-destruction process is two or more, 
then it is highly unlikely for the reverse process to occur. 
The low rate for the reverse process leads to the absence 
of a bottleneck. 

On the other hand, phonon-mediated recombination 
is inherently susceptible to a bottleneck because only 
one particle (the phonon) is created when the pair is de- 
stroyed. Following pulsed injection, the bottleneck sets 
in when the pair destruction and creation rates become 
equal, or Pn^^ = 2jpcNss, where nss and Ngs are the 
quasiequilibrium steady-state densities of quasiparticles 
and 2A phonons, respectively. The size of nss relative to 
the initial density no is a measure of the strength of the 
bottleneck. This ratio can be calculated using conserva- 
tion of energy, which implies that 2nss + Nss — 2no. The 
result is 

no 2 V A 

where A = 2"fpc/ f3nQ. The strong bottleneck regime cor- 
responds to A ^ 1 , in which case the steady state density 
nearly equals the initial density, nss /no ~ 1 — 1/A. On 
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the other hand, if A 4; 1 the bottleneck is practically 
nonexistent and the steady state density is very small, 

To determine whether the cuprates are in the strong or 
weak bottleneck regime, we must estimate the magnitude 
of A. To do so we first recognize that A must exceed a 
Amin, obtained by setting no equal to iV(0)A, the maxi- 
mum quasiparticle density in the superconducting state. 
We can evaluate Xmm in the dirty s-wave case, as Ka- 
plan et al. [sJl have calculated as well as (3 in this 
regime. The characteristic phonon decay rate is given by 
7pc = 47r^iV(0)a^A//i7V, where N is the number density 
of ions. Given this phonon decay rate, we find that 



27r 7V(0) 
Y NF{2A)' 



(12) 



The simplicity of this result is a consequence of detailed 
balance in the dirty- limit. Eq. 12 shows that the strength 
of the bottleneck is determined by the ratio of the elec- 
tron density of states at the Fermi level to the phonon 
density of states at 2 A. It is clear that the high-Tj, 
cuprates can be in the weak bottleneck regime, where 
A is very small, in constrast with conventional low-Tc 
superconductors. In the high- Tc materials the gap en- 
ergy is large, and 2A coincides with large phonon density 
of states peaks associated with optic phonons. Further- 
more, the density of states at the Fermi level is small com- 
pared with conventional metals. Both factors favor small 
A. In contrast, in conventional superconductors phonons 
with energy 2A lie in the low-energy tail of the acous- 
tic spectrum. Therefore the density of resonant phonon 
states is small. This fact, together with the large elec- 
tronic density of states, place conventional superconduc- 
tors in the strong bottleneck regime. 



VI. SUMMARY 

The most significant observation reported in this paper 
is the strong dependence of the decay rate of the tran- 
sient reflectivity on both temperature and pump laser 
intensity. At low T the decay rate is a linear function of 
I. At low / the decay rate decreases exponentially with 



decreasing T. At the lowest values of T and / probed 
by our experiments, the lifetime of the photoexited state 
becomes extremely long, « 600 ps. There is no indica- 
tion that the lifetime would not continue to increase if 
the experiment could access lower values of T and /. 

In our analysis of the data, we first addressed the na- 
ture of the excitations probed by measuring AR. We 
presented arguments that AR is proportional to the den- 
sity of high-energy or antinodal quasiparticles created 
by the laser excitation. We explained theoretically why 
such particles, although not the lowest energy excita- 
tions, could be highly metastable at low temperature and 
density. The metastability is a consequence of the limited 
phase space for decay of high-energy quasiparticles to the 
nodal regions of the Brillouin zone. As the phase space 
restrictions rely on momentum conservation, the degree 
of metastability may be expected to be highly sensitive 
to disorder. 

Having understood the stability of an isolated quasi- 
particle, we sought a quantitative understanding of the 
decay process that occurs when a pair of quasiparticles 
interact. We estimated the magnitude of the recombi- 
nation coefficient, /?, from both the T and / dependence 
of the decay rate. We emphasized that /? is a direct 
measure of the strength of quasiparticle interactions. In 
addressing the theoretical approaches to the recombina- 
tion rate, we identified h/m* as the natural unit for (3 in 
a two-dimensional superconductor. We pointed out that 
while recombination is assumed to be phonon-mediated 
in conventional superconductors, electron-electron medi- 
ated processes must be considered as well in the cuprates. 
In the context of phonon-mediated processes, we dis- 
cussed the need to extend the work of Kaplan et al. [sj 
to the 2D, clean limit that is applicable to the cuprates. 
In the electron-electron mediated case, more work is 
needed to understand the underlying recombination pro- 
cess. As a motivating factor, we pointed out the remark- 
able coincidence between the normal state scattering rate 
and the superconducting state recombination rate. We 
showed that both processes have the same dependence 
on the quasiparticle density, which can be described by 
an interaction coefhcient equal to H/Trm*. 

This work was supported by NSF-9870258, DOE-DE- 
AC03-76SF00098, CIAR and NSERC. 
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